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We present new evolution sequences for very low mass stars, brown dwarfs 

^ ■ and giant planets and use them to explore a variety of influences on the evolution 

of these objects. While the predicted adiabatic evolution of luminosity with time 

is very similar to results of previous work, the remaining disagreements reveal the 

magnitude of current uncertainty in brown dwarf evolution theory. We discuss the 

sources of those differences and argue for the importance of the surface boundary 

• condition provided by atmosphere models including clouds. 

The L- to T-type ultracool dwarf transition can be accommodated within the 
00 ■ I " ^ i i 

O ■ lAckerman fc Marleyl (120011 ) cloud model by varying the cloud sedimentation pa- 

00 ' 
O 



rameter. We develop a simple model for the evolution across the L/T transition. 
By combining the evolution calculation and our atmosphere models, we generate 
colors and magnitudes of synthetic populations of ultracool dwarfs in the field 
and in galactic clusters. We focus on near infrared color-magnitude diagrams 
(CMDs) and on the nature of the "second parameter" that is responsible for the 
scatter of colors along the T e g sequence. Instead of a single second parameter we 
find that variations in metallicity and cloud parameters, unresolved binaries and 
possibly a relatively young population all play a role in defining the spread of 
brown dwarfs along the cooling sequence. We also find that the transition from 
cloudy L dwarfs to cloudless T dwarfs slows down the evolution and causes a pile 
up of substellar objects in the transition region, in contradiction with previous 
studies. The same model is applied to the Pleiades brown dwarf sequence with 
less success, however. Taken at face value, the present Pleiades data suggest that 
the L/T transition occurs at lower T e g for lower gravity objects, such as those 
found in young galactic clusters. The simulated populations of brown dwarfs 
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also reveal that the phase of deuterium burning produces a distinctive feature in 
CMDs that should be detectable in ~ 50-100 Myr old clusters. 

Subject headings: stars: low mass, brown dwarfs — stars: evolution — stars: 
atmospheres 



Introduction 



There are now approximately 450 L dwarfs and 100 T dwarfs known (see iKirkpatrick 



(120051 ) for a review of these spectral classes). They span effective temperatures from about 
2400 to 700 K and exhibit a range of gravities, metallicities, and atmospheric condensate con- 
tents. After more than a decade of intense study, the modeling of the complex atmospheres 
and synthetic spectra of brown dwarfs has reached a rather high degree of so phistication, 



including the chemis try of a very large number of gas and condensate species f Allard et al. 



2001 



2008 



2002 



2005) , and particulate c loud models (lAUard et al.l l200ll ; lAckerman &: Marleyl l200ll ; iTsuji 



.Lodders 2002), i ncrea singly complete molecular opacity databases fjFreedman et al. 
Sharp fc Burrows 2007 ). extreme resonance line broadening feurrows et al. Il2000l : lAUard et al. 



Helling et al.l 120081 ). While a few conspicuous problems remain, the synthetic spectra 
and colors reproduce the observations fairly well and the determination of the basic astro- 
physical properties of brown dwarfs has begun in earnest. The full astrophysical benefit 
of synthetic spectra and colors is obtained when atmosphere calculations are coupled with 
evolution models that provide the surface atmospheric parameters (T e s,g) as a function of 
mass and age. The time evolution of the spectrum and colors, as well as absolute fluxes can 
be computed and directly compared with observations to estimate astrophysical parameters 
that are not easily amenable to direct observation. 

To enable such comparisons using our own model atmosphere effort and to pursue a more 
complete analysis of spectroscopic and photometric data, we have developed a code to com- 
pute evolution sequences of low mass stars, brown dwarfs and giant planets. These evolution 
seque nces have already been applied extensiv e ly to the analysis of brow n dwarf observa- 
tions jRoellig et al.l l2004j : ISaumon et all bood . I2OO7I : iLeggett et aD l2Q07al fbL ICushing et al 



20081 ) but have not been discussed in any detail. Here, we describe the input physics and 
assumptions of our particular approach as well as the unique aspects of our atmospheric 
boundary condition. The evolution of isolated brown dwarfs — a relatively simple case of 
stellar evolution — has been studied extensively for the past 20 years and is well understood 
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is similar to the more recent work and the result are very similar to previous work. A de- 
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tailed comparison with other published calculations of the evolution of brown dwarfs reveals 
small differences that we quantify and, to the extent possible, attribute to the different as- 
sumptions and approximations in each model. We highlight the application of the evolution 
sequences to the calculation of synthetic near infrared color-magnitude diagrams (CMDs) 
to explore the nature of the "second parameter" responsible for the spread of observed ob- 
jects around the main trends along the L and T spectral sequences. We develop a simple 
parametric model for the L/T transition to reproduce with good success the CMD of field 
brown dwarfs. We apply the same model to synthesize the population of the Pleiades clus- 
ter (110 Myr) and compare with the latest deep survey data for this galactic cluster. We 
discuss several potentially observable features in near-infrared CMDs that would illuminate 
the evolution of brown dwarfs as well as the nature of L/T transition. 



2. Evolution model: Assumptions and physical inputs 

The evolution model assumes adiabatic cooling of spherical, hydrostatic, non-magnetic, 
non-rotating brown dwarfs. The adiabatic assumption is valid for dense fully convective 
structures and objects with masses ranging from 0.3 M & down to Saturn's mass can be 
modeled. In this section we consider key physical aspects of the evolution model. 



2.1. Electron Conduction 



The effects of electro n conduction on the cooling of old brown dwarfs, pointed out 
by IChabrier et al.l (j2000al ). are not included in our adiabatic calculation. Conduction be- 
come s significant only for t he more massive brown dwarfs and at ages greater than about 2 
Gyr (IChabrier et al.l l2000al ) . Because the structure of a brown dwarf model is determined 
primarily by the pressure of the degenerate electrons but its heat content is given by the 
temperature of the ions, a cooling brown dwarf goes through very nearly the same sequence 
of L,T e s, R, but at later times when heat transport by conduction is included. We will see 
in §4 that neglecting conduction produces a relatively modest error when considering other 
sources of uncertainty. 
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2.2. Thermal Structure and evolution 



The adiabatic evolution is obtained by solving the mass equation 



dm 
dr 



the equation of hydrostatic equilibrium 

dP 

dr 

and the equation of conservation of energy 



47rpr^ 



pGm(r) 



M 



e nuc dm I dt 



■)■ 



M 



TdS dm, 



(2) 



(3) 



where the rate of nuclear energy generation is e nuc (p, T) and the equation of state (EOS) 
along an adiabat is P(p, S) and T(p, S). In these equations, p, P, T and S are the density, 
the pressure, the temperature, and the entropy of the gas, respectively, L = 4iiR 2 aT^ s is 
the luminosity, and m{r) is the mass interior to radius r. The equation of conservation of 
energy controls the evolution time scale as it gives the time step dt between two successive 
adiabatic structures that differ in entropy by an amount dS. 

We use the hydrogen and helium EOSs of Saumon, Chabrier & Van Horn (1995, hereafter 
SCVH), which was developed specifically for this type of application. This EOS is common 
to nearly all other evolution and structure calc ulations of low mass stars, b r own dwarfs, ex- 



(2000a 


b); 


Baraffe et al. 


(2002 





. Appendix A reports on typo- 
graphical errors that were recently found in SCVH. The contribution of the metals to the 
EOS, which results in a slightly smaller radius for given mass and entropy, is neglected. 



2.3. Nuclear Energy Generation 



Nuclear energy generation in low mass stars and brown dw arfs is quite simple and is 
reduced to one branch of the pp chain (IBurrows fc Liebertlll993l ): 



and 



p + p — > d + e + + v e 
p + d -> 3 He + 7. 



(4) 
(5) 
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We use the nuclear reaction cross-sections from the NACRE data base (lAngulo et al. 



19991 ) and apply the screening corrections developed by Chabrier which include both ionic 



and electronic screening (Saumon et al. 1996; Chabrier & Baraffe 1997; Chabrier, priv. 
comm.). We assume that the compositional changes due to nuclear reactions are homog- 
enized by convection throughout the brown dwarf on a time scale much shorter than the 
nuclear burning time scale (see Chabrier fc Baraffe 1997, however) . Lithi um burning is a 
useful age/mass diagnostic of brown dwarfs (IRebolo et al.lll992l ; lBasrill2000l ) but contributes 
negligibly to nuclear energy production and is not included in our code. 



2.4. Initial State 

The initial state for the evolution assumes an extended, high entropy, spherical config- 
uration defined by T e g > 4000 K (higher for larger mass). Such a n idealized init i al con dition 



bears little relation with the actual formation of brown dwarfs. iBaraffe et al.l (120021 ) have 
shown that the evolution of brown dwarfs is quite sensitive to the choice of initial state for 
ages under ~ 1 Myr. For this reason, we present only the evolution for later times. The 
early evolution of our models is further limited by the high-temperature limit of the atmo- 
sphere grid (see below). For the initial interior composition, we adopt Y = 0.28, Z = 0, and 
D/H= 2 x 10~ 5 by number (2.88 x 10 -5 by mass), and we set the abundance of 3 He to zero. 



2.5. Surface Boundary Condition 



The most significant difference between the BD evolution calculations published over the 
past decade is in the treatment of the surface boundary condition of the model, which con- 
nects the surface properties (T e g,g) to the interior model (M, L, R, age). For our adiabatic 
(i.e. constant entropy) interior models, applying the surface boundary condition consists of 
matching the interior entropy to that of the convective bottom of the atmosphere model. 
The importance of using realistic, non-gray atmosphere models for the surface boundary 
condition of very- low mass stars and brown dwarfs has been emphasized for some time 
( jSaumon et al.lll994t IChabrier fc Baraffeill997l ). We use our grids of atmosphere models to 
obtain the surface boundary condition which is expressed as the entropy at the bottom of 
the atmosphere as a function of the atmospheric parameters, S(T e g,g, [M/H], / se d), where 
f se d is a cloud model parameter (see below). This makes the evolution fully self-consistent 
with the modeled spectra. 

Our atmosphere models have been described previously (IMcKay et al. Il989l : iMarley et al 
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An updated, detailed description will be the subject of an upcoming publication. Briefly, the 
model computes radiative-convective equilibrium atmospheric temperature-pre ssure struc- 



tures . Our thermal radiative transfer follows the source function technique of iToon et al. 



( 119891 ) allowing inclusion of arbitrary Mie scattering particles in the opaci ty of each layer. 



Our o pacity database, accounting for all important absorbers is described in lFreedman et al. 
( 120081 ). Our chemical equilibrium grid of molecular, atomic and ionic abu ndances as a func- 



tion of tem p erature, p ressu r e, an d metallicity is bas e d on the work of iFeglev fc Lodders 
dl994j . Il996l ): lLoddersI Jl999l . I2OO2I) : lLodders fc Fedevl J2OO2I ) and lLodders fc Fedevl J2006l)~ 
We u se the elemental abundances of lLodderd (120031 ) . Our baseline cloud model (jAckerman &: Marlev 
20011 ) parametrizes the efficiency of sedimentation of cloud particles through an efficiency 
factor, / se d. Large values of / se d correspond to rapid particle growth and large mean particle 
sizes. In this case condensates quickly fall out of the atmosphere, leading to physically and 
optically thin clouds. When / sec j is small, particles grow more slowly and the atmospheric 
condensate load is larger and clouds thicker. Our cloud model is fully coupled with the 
radiative transfer and the (P, T) structure of the model during the calculation of a model so 
that they are fully consistent when convergence is obtained. The mod el has been very suc- 



cessor 



2001 



when applied to the atmospheres of brown dwarfs and Jupiter (jAckerman fc Marlev 
Knapp et all 12004 : iGolimowski et al.l I2004J ). In particular Marley et al. (2004) and 



Cushing et al.l (120081 ) show generally excellent fits between our model spectra and observa- 
tions of cloudy L dwarfs. The near infrared colors of brown dwarfs are quite sensitive to the 
choice of / se d, a point we will return to in the discussion of color-magnitude diagrams (§4). 

Our grid of atmosphere models covers 500 < T c g < 2400 K and 3.5 < log g < 5.5, which 
does not provide a boundary condition for the late evolution of low-mass objects and the very 
early evolution of the more massive ones. To obtain a boundary condition for T e g < 500 K, we 
define Tm as the temperatu re at P = 10 bar that gives the same entropy S as the atmosphere 
model (ILunine et al.lll989l ). We smoothly interpolate Ti {T e s,g) quadratically to T w = 
at T e ff = at constant gravity. The interpolated values of T 10 are then converted back to 
an equivalent entropy S(T W ,P = 10 bar) for use as the surface boundary condition. For 
T e ff > 2400 K and gravities outside the range of the atmosphere grid, linear extrapolation in 
T e ff and log g are used. This gives well-behaved but inaccurate results for the early, warm 
phase of the evolution and for very-low mass objects (M < 0.002 M ). 

In the following we present two evolution sequences, one based on a cloudless atmosphere 
grid and one using a cloudy grid with / se d = 2, both with solar metallicity. The two surface 
boundary conditions (expressed as T10) are shown in Fig. [1] for selected gravities. Clouds 
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play a minor role in high-T e ff atmospheres and both sequences converge to the same boundary 
condition. The cloudy surface boundary condition is not as smooth as the cloudless case 
because of occasional numerical difficulties in the cloudy atmosphere calculation. Cloudy 
atmospheres have a globally larger opacity than clear models, which results in higher entropy 
in the convection zone (i.e., a hotter interior structure, or higher Tio, Fig. [1]). The differences 
between the cloudy and cloudless surface boundary conditions are significant. Calculations 
of the evolution of cloudy models with other values of / se d are possible and desirable but 
require an increasing number of model atmospheres to define the surface boundary condition. 
The boundary conditions from restricted grids of models for / se( j = 1-4 shows less than half 
of the difference in T 10 that is seen between / sc< j = 2 and cloudless models, and even less in 
some regimes. To a good approximation, the / se d = 2 evolution sequence can be used for all 
cloudy cases. 



3. Evolution sequences 



3.1. Cloudless and cloudy surface boundary conditions 

The evolution of the luminosity is shown in Fig. [2] for cloudless models and Fig. [3] 
for the cloudy models for objects ranging from 0.005 to 0.08 M & . Both show the well- 
known features of the evolution of very low mass star s, brown dwarfs and giant planets 
(ID'Antona fc Mazzitellil Il985l ; iBurrows et al.lll989l . 119971 ). with an initial deuterium burning 
phase (M > 0.015 M Q ), followed by contraction and eventual settling on the main sequence 
for the higher mass models (M > 0.075 M ^). The effect of atmospheric opacity on the 



cooling of brown dwarfs is also well known (Burrows et al.lll993l; IBurrows fc Liebertl Il993 



Chabrier fc BaraffelboOO ; Ichabrier et aliboOOal ) and a comparison of the cloudy and cloudless 
sequences shows the expected trends. The cooling rate is primarily controlled by the slope 
of the T 10 (T eff ) relation (F ig. [Tj and for T eff > 1 000 K, the cloudy models are cooler and less 
luminous for a given age (jChabrier et al.ll2000al ). We find that at T eS = 1400 K, a cloudless 
brown dwarf is ~ 35% older than one with a / se d = 2 cloudy atmosphere, a value that is 
nearly independent of mass from 0.02 to 0.06 M Q . 

For these evolution sequences with solar metallicity atmospheres, the minimum hydrogen- 
burning mass, defined by objects that reach stable H-burning at an age of 10 Gyr is 0.075 M Q , 
where T e fj = 1910 K for the cloudless sequence. The main sequence of cloudy models ends 
at a lower mass of 0.070 M Q and T e g- = 1550 K. The latter falls within the range where 
clouds have a maximal effect on the atmosphere and is therefore a more realistic value of the 
hydrogen minimum burning mass of solar metallicity stars than the cloudless case. Defining 
the minimum deuterium-burning mass as objects that burn 90% of their initial deuterium in 
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10 Gyr, we find 13.1 Mj upiter and 12.4 Mj upiter for the cloudless and cloudy cases, respectively. 

For comparisons with observables, a representation of the evolution in terms of T e g and 
gravity is more transparent (Figure H]). The range of T e ff shown covers all spectral types later 
than L0, and includes the bottom of the main sequence. The figure also shows isochrones, 
and lines of constant radius and luminosity. The knowledge of any two quantities among 
T e ff , g, R, M, L and age allows a determination of all the others. The most notable feature is 
a bump in all isochrones for masses between 0.01 and 0.02 M which signals the deuterium 
burning phase. This figure shows a number of interesting features of the evolution. For 
example, an object with a known L (red curves) has an upper limit to its T e g. The figure 
also shows that there is a maximum gravity for any brown dwarf in a sequence of models. 
This arises from the mass-radius relation of old, degenerate brown dwarfs which is primarily 
determined by the equation of state. For the cloudless sequence of Fig. HJ the upper bound 
on the gravity is logg = 5.465, which occurs for a 0.068 M & brown dwarf at T e g = 1160K. 
In the cloudy sequence, the peak gravity occurs at the same mass but logg = 5.366 and 
T e fj = 1380 K. There is a similar weak dependence of the maximum gravity on the metallicity, 
which increases by 0.023 dex for a decrease of 0.3 dex in [M/H] (for cloudless sequences). 
More generally, there is an upper (T eff ,logg) envelope to the evolution of brown dwarfs that 
in practice is defined by the main sequence and the oldest objects (~ 10 Gyr for the disk 
population). This envelope puts a firm upper limit on determinations of the gravity of brown 
dwarfs (for a given T e g) such as those obtained by fitting spectra. Higher values of the gravity 
are unphysical. 

The (Tgff , log g) evolution of the cloudy sequence is shown in Fig. [5] where the cloudless 
cooling tracks are also shown for comparison. The effects of clouds on the surface boundary 
condition (Fig. [TJ essentially vanish at high and low T e g , and consequently, the cooling tracks 
of the cloudy models overlap the cloudless tracks in both limits. At intermediate T e g, where 
the cloud effects are maximal, the gravity of an object of a given mass is up to 0.08 dex lower 
in the cloudy models than for the cloudless case. For given values of (T e g,\ogg), differences in 
log L can range for to 0.1 dex, from to 6% in radius, and to 60% in age. These differences 
are modest but will become significant as brown dwarf studies become increasingly detailed 
and precise. 



3.2. Comparison with other calculations 

We validate our calculation by comparing with published results that are widely cited 
in the brown dwarf literature. This also offers a measure of the relative uncertainties that 
remain in modeling the evolution of brown dwarfs and, to some extent, to identify the model 
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assumptions responsible for t hose differences. Fig ures [2] and [3] compare our cloudless sequence 
with the COND03 mo dels of IBaraffe et al.l (120031) and our cloudy f KP ^ = 2 sequence with the 



DUSTY00 models of IChabrier et al.l (j2000a[ ) and lBaraffe et al.l (J2002|), respectively. We note 
that the COND03 and DUSTY00 models and ours use the same H/He EOS, the same interior 
composition and that there is no significant difference in our respective nuclear reactions 
rates (including screening factors). Except for the modest role of electron conduction, the 
differences can only arise from the initial boundary condition (for ages under a few Myr) 
and the surface boundary condition provided by the atmosphere at later times. 

Generally, the COND03 and DUSTY00 luminosities at a given age are slightly higher 
than ours. In both cases, however, the agreement is excellent for all masses at all ages above 
a few Myr, with differences below 0.05 dex for the cloudless case and below 0.09 dex (and 
usually much less) for the cloudy sequence. At young ages, the differen ces can be significan t 
and reflect our different choice for the initial conditions for the models (IBaraffe et al.ll2002l ). 
For massive brown dwarfs at late times, electron conduction becomes the dominant energy 
transport mechanism in the core (§2.1). At 10 Gyr of age, the luminosity of the COND03 
models < 0.06 M is systematically higher than ours, an effect that decreases for lower 
masses. This is the expected behavior for models that include electron conduction. For a 
10 Gyr old brown dwa rfs of 0.06 M Q , we find a luminosity that is 0.104 dex l ower than that 



of B araffe et al.l (120031 ) , which is in perfect agreement with the value given in IChabrier et al. 



( 12000al ) for the change in luminosity brought by introducing conduction in the evolution. 
Thus, the differences at late times are primarily due to our neglect of conduction and those 
differences are rather small. 



A comparison of our cloudless models with the evolution sequence of iBurrows et al. 



( 119971 ) (hereafter, B97) is shown in Fig. [6j The agreement is not as good as we found 
with the COND03 models (Fig. [2]). One of the reasons is that the B97 models use a 
rather low value for the helium abundance of Y = 0.25, while the present mo dels and 
the C OND03 models use Y = 0.28. The proto-solar value is Y = 0.274 ± 0.12 jLodders 
20031 ). To better compare with B97, we have computed a short sequence of models with 
Y = 0.25. This systematically shifts the luminosity curves downward and reduces the 
differences but not enough to reach the level of agreement seen in Fig. [2J The B97 models 
also use the same EOS and nuclear reaction rates as in our calculation. The screening 



correction (IGraboske et al.lll973l ) is more crude but this only affects models at the edge of 
the main sequence. The primary reason for the differences, and for the structure seen in the 
B97 L(t) curves at logL/L ~ —4 that is not apparent in the COND03 or in our cloudless 
sequence is the surface bou ndary condition. Th e model atmospheres used in B97 are the 



non-gray cloudless models o f Mar 



at the higher T e g (jSaumon et al. 



ev et al.l (119961 ) . complemented with gray cloudless models 
19961 ). Atmosphere models of brown dwarfs have become 
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increasingly more realistic over the past decad e. Particularly significant changes in cloudless 
models since the work of iMarley et al. are 1) the recognition of th e role of the K I 

and Na I reso nance doublets at optical wayelengt hs (IBurrows et al.ll2000l ). 2) improved 



CIA opacity (Borvsow et al. 2001 ; Borysowl I2OO2I ) . 3) a new line list for TiO ( Allard et al. 



2000 ; Freedman et al. 20081 ). 4) improved modeling of the condensation che mistry (Lodders 



2002]), 5) new sources of molecular opacity such as CrH, FeH, VO and PH 3 ( jFreedman et al 



( 120081 ) and refe rences therein), and 6 ) expanded molecular opacity lines list - notably for 
CH 4 and NH 3 ( jFreedman et al.l 120081 ) . While the evolution is not sensitive to the details 
of the atmospheric opacity, the global opacity of the atmosphere largely accounts for the 
differences between the various evolution calculations for brown dwarfs. The above factors 
and cloud formation as well, all play a significant role in increasing the a tmospheric opacity. 



each i n a different range of T e fj and gravity. The analytic evolution of IBurrows fc Liebert 



( 119931 ) show that the luminosity increases with the Rosseland mean atmospheric opacity as 
relation confirmed by our non-solar cloudless evolution sequences. This, along 
with the choice of Y — 0.25 explains most of the differences between our calculation and 
that of B97. 



4. Synthetic color-magnitude diagrams 



The ayailab ility of parallaxes for brown dwarfs (IDahn et al.l |2002| ; iTinney et al.l 12003 



Vrba et al.l 120041 ) provides absolute magnitudes and spectral fluxes that can be compared 
with a combination of evolution models and synthetic spectra and colors. The near infrared 
CMD of field brown dwarfs (Fig. [7]) shows a number of remarkable features. The L dwarf se- 
quence naturally extends the M dwarf main se quence, where evidence for condensates in the 
atmosphere first appears (IJones fc Tsujilll997l ) suggesting that clouds play a continuous role 
from late M dwarfs through the L spectral sequence. The photometric sequence is increas- 
ingly broader with later spectral types, up to the latest L dwarfs, which su ggests the presence 
of a '' se cond paramet e r" be sides T pff that affects the co lors of L dwarfs. iTsuji fc Nakajima 



(120031 ). iKnapp et al.l (120041 ) and IBurrows et al.l (120061 ) have suggested that variations in 



gravity (i.e. a spread in mass and age) account for the broadening. Other factors such as 
metallicity, unresolved binaries, and a range of cloud properties may also contribute to the 
scatter. Between the latest L dwarfs and the mid-T dwarfs, there is a dramatic shift of 
J — K to the blue, long recognized as due to the lack of cloud opacity in late T dwarfs and 
the appearance of CH 4 in the K band. The late T dwarf sequence is also fairly broad, pro- 
viding further evidence of the importance of additional physical parameters in brown dwarf 
atmospheres. An important revelation enabled by C MDs is that mid-T dwarfs are brighter 



in J (but not in H or K) than the latest L dwarfs (IDahn et al.l 120021 ; ITinney et al.l 12003 
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Vrba et al. 



200 4]), even though the bolometric luminosity decreases steadily with later spec- 



tral types (IGolimowski et al.l 120041 ) . This was confirmed by the discovery of binary brown 



dwarfs w ith (coeval) components spanning the L/T transition and showing the J band flux 
reversal (IBurgasser et al.ll2006l ; iLiu et al.ll2006l ). Even more intriguing, the binary fraction of 
early T dwarfs is twice as large as for other dwarfs and their spectra and colo rs are largely the 



resul t of the combined light of components of earlier and later spectral types (IBurgasser et al. 



20061 ; ILiu et al.l 120061 ) . Together, these facts bear on the nature of the transition from L to 
T spectral types which is believed to be caused by a rather dramatic change in the cloud 
properties (such as the fraction of cloud c overage, particle s ize, and the cloud's ver tical ex- 
tent) over a small range of T e g- of ~ 200 K (I Vrba et al.l 12004 ; iGolimowski et al.ll2004l ). There 
is as yet no physical modeling of such a process. 

Detailed comparisons of our self-consistent ev olution and spectra with sp ectroscopic data 
of a few late T dwarfs have been rather successful (jSaumon et al.ll2006l . 120071 ) , demonstrating 
that at least for the coolest brown dwarfs known the models, while admittedly imperfect, 
are fairly realistic. Fits of the entire spectral energy distribution of brown dwarfs with 
earlier spectra types are more challenging due to the cloud opaci ty but the models h ave 
reached a level of realism where meaningful results can be obtained (jCushing et al.ll2008l ). A 
comparison of our model colors in the mi d- and near-infrared w ith a sample of ~ 50 brown 
dwarfs shows generally good agreement (ILeggett et al.l l2007bl ). The near infrared CMDs 
of our model colors are shown in Fig. [8] which displays the trends with T e s, gravity and 
the cloud sedimentation parameter. The / se d = 1 sequence generally follows the observed L 
dwarf sequence and steadily becomes redder with decreasing T e fj0 In models with decreasing 
cloud thickness (increasing / se d), the J — K color is the same as for / sec j = 1 at the high-T e g 
end where the cloud opacity is always modest but shows a turnover to the blue that occurs at 
gradually higher T eS as the thinner cloud deck sinks below the average photosphere. At low 
T eff , the colors of / se d = 4 models join those of cloudless models (/ se d — > oo) below ~ 800 K. 
The cloudless models do not show the redward trend of the L sequence but grow steadily 
bluer in near-infrared colors, due to the onset of H 2 collision-induced absorption and of CH 4 
bands. 

In a very qualitative sense, the cloudy models show the turnover in J — K seen in the 
L/T transition objects but the turnover is much too gradual. Furthermore, the cloudy mod- 
els do not show the brighte ning in the J band. A com parison with Fig. [7] indicates that 
within the framework of the lAckerman fc Marleyl (120011 ) cloud model, the L to T transition 
can be modeled by increasing / se d during the cooling of a brown dwarf. This idea is sup- 



1 We note however that f se d = 1 models prod uce relatively flat-topped JHK band peaks, a feature not 
seen in spectroscopic data ((Stephens et al.ll2008l ). 
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ported by the detail ed analysis of the spectral energy distribution of several brown dwarfs 
( ICushing et al.ll2008l ) . While it is plausible that the transi tion could be caused by an i ncrease 
in sedimentation efficiency with T e g at the transition, the lAckerman &: Marleyl (120011 ) model 
does not provide a physical justification for this process. The transition could be the conse- 
quence of other cloud processes that are not captured in the Ackerman Sz Marley model. The 
rather simplistic cloud models developed so far are not a dequate to understan d the cause of 
the transition and the complex dust formation model of iHeUing et al.l (120081 ) has not been 
applied to modeling the L/T transition. Despite these caveats, the one-parameter Ackerman 
&: Marley cloud mod el reproduces observed spectra of brown dwarfs of all spectral types 
(ICushing et al.l 120081 ) rather well and provides a useful basis for further analysis of data. 
Cloud models with more detailed physics, while weakly constrained by present observations, 
are very much needed however. 

Color-magnitude diagrams have historically been one of the most powerful tools to 
understand stellar evolution. The difficulty in discovering and observing the intrinsically 
faint brown dwarfs has limited the the application CMDs in this brand of stellar astrophysics. 
Yet, a comparison of modeled magnitudes and colors with a large volume-limited sample of 
brown dwarfs with measured parallaxes could enable the determination of the cooling rate as 
a function of T eff , of the T e g range of the transition from cloudy to cloudless atmospheres, the 
mapping of the variation of cloud parameters with T e g, and reveal the distributions of cloud 
parameters and gravities at a given T c g. Such a sample is not yet available. Furthermore, 
the models are not quite reliable enough to carry out such a program. On the other hand, 
the colors of synthetic populations of field brown dwarfs is a tool to explore the nature of 
the "second parameter" along the L dwarf sequence that complements the clues gathered 
from brown dwarf binaries. In this section, we explore the CMDs of simulated brown dwarf 
populations based on our evolution sequences and model atmospheres, compute a hybrid 
evolution sequence that models the L/T transition and make qualitative comparisons with 
data in the M K vs J — K CMD to identify trends under various assumptions. We first 
consider the population of field brown dwarfs, followed by a discussion of brown dwarfs in 
the Pleiades. 



4.1. Disk brown dwarfs 

We simulate a volume-limited population of field brown dwarfs by generating a random 
set of points in the mass-age parameter space with the mass chosen between 0.006 and 
0.1 M , assuming a power law initial mass function (IMF) 
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where dN is the space density of brown dwarfs in mass interval dM. We choose a — 1, 
which is represe ntative of the mas s function of field brown dwarfs and consistent with the 
determination of lAllen et al.l (120051 ). We adopt a constant star formation rate (SFR) over the 
age of the Milky Way (0 - 10 Gyr). The normalization of the number density iV is irrelevant 
for our purposes. We arbitrarily normalize all simulations to a fixed number of objects in the 
mass range 0.075-0.08 M®. Because all brown dwarfs with known parallax are nearby, we 
neglect extinction and reddening. For simplicity, our synthetic population consists of single 
dwarfs only. The effects of binaries and of different assumptions for the SFR and the IMF 
are discussed separately in §4.2.1. 

We convert our sample of (M, age) points to (T eff ,logg) values with our evolution 
sequences. This is shown in Fig. M where we have used the cloudy evolution sequence with 
f se d = 2. Objects with T e g above the upper limit of our surface boundary condition (2400 K) 
have been culled from the sample. The distribution can be compared with the cooling 
tracks and isochrones for this sequence (Fig. [SJ). The resulting (T cS , logg) distribution is 
rather remarkable. As can be expected for a population with an average age of 5 Gyr, the 
vast majority of objects are concentrated near the outer envelope of the evolution defined 
by the 10 Gyr isochrone. The upper boundary, between logg = 5.3 and 5.4, contains a 
narrow distribution formed by the low-mass end of the main sequence and the most massive 
brown dwarfs (M > 0.06 M ). The vast majority of objects are found at T c s < 1000 K; 
a consequence of the rising IMF at low masses and of their rapid cooling. Changing the 
IMF power-law index a alters the relative abundance of points along the low-T e g edge of the 
distribution but the general appearance of this figure is unaffected: A thin distribution at 
high gravities and a broader, more populous distribution at low-T e ff. Most of these objects 
are very dim, however, as logL/L < —6 for T e s < 600 K. Figure [9] illustrates that warm, 
low gravity brown dwarfs would be quite rare in the field under the assumption of a constant 
SFR rate over the age of the Galaxy. In contrast, t he kinematics of L and T dw arfs indicates 
that the population may be as young as 0.5-4 Gyr (jZapatero Osorio et al.ll2007l ). a possibility 
we consider in 54.2.1. 



The pi o neerin g wor k of IChabriei (120021) and the extensive studies of iBurgasserl (120041 ). 
Allen et al.l J2003h and kllen et al.l~(l2005h show how variations in the IMF and SFR af- 
fect one-dimensional distribution functions (such as the luminosity function) of field brown 
dwarfs. They offer a clear and detailed discussion of the features of the very-low mass stellar 
and substellar luminosity function in terms of the characteristics of their evolution. These 
studies focus primaril y on a determination of the substellar mass function and the mass 
budget of the Galaxy (jChabrierl 1200 ll . 120021 ) or on the luminosity function which collapses 
all objects from a two-dimensional distribution - for example, age and mass - into a one- 
dimensional distribution, binned in terms of T e R , L, or spectral type. Our primary concern is 
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to explore the two-dimensional distribution and its potential to reveal additional features of 
the evolution of brown dwarfs by constructing color-magnitude diagrams. Figures (fTUl) show 
our synthetic field distribution of brown dwarfs in the (Mr, J — K) CMD for different values 
of the cloud condensation parameters / se d- All four figures are based on the same cloudy 
evolution sequence with / se d = 2 and the (T e g,g) distribution of Fig. [9j but each uses colors 
computed from atmosphere models with different values of / se d- An evolution sequence with 
a single value of / se d should be adequate for most purposes (see §2.5). Finally, Fig. [TUfe 
shows the cloudless case, computed with the cloudless evolution sequence and cloudless at- 
mosphere models. The figures also shows the photometry of M,L, and T dwarfs with known 
parallaxes, excluding known binaries with blended photometryo We note that since the set 
of dwarfs with measured parallaxes is not in any sense a volume-limited sample and is sub- 
ject to unquantifiable biases, the comparisons must be regarded with caution. For the same 
reasons a detailed, quantitative comparison of the modeled and observed two-dimensional 
distributions would be premature and our study is mostly qualitative in nature. 

The most striking feature in these diagrams is that the color-magnitude distribution 
of the sample is rather different from that of the simulation, the former being more or less 



(Allen et al. 


2005; 


Bureasser 


2004) 



This shows that the sample of brown dwarfs with measured parallaxes is more similar to a 



magnitude limited sample than the volume limited sample we simulated (jChabrierll2002l ) . We 
now turn to a comparison of our synthetic CMDs with the observed features, starting with 
the high-T e fj end of the sequence and proceeding to the cooler brown dwarfs. Simulations 
under different assumptions are discussed in §4.2.1. 

At the upper limit of our calculation (T eff = 2400 K), the synthetic sequence joins the 
late M sequence nicely. The cloudless and cloudy sequences meet at T eff ~ 2400 K as clouds 
play a minimal role at higher T e g (Fig. [H]). In the synthetic population, the J — K color 
of the objects with 10.5 < Mk < 11-5 becomes too blue by ~ 0.3, a problem that has 
plagued low temperature atmosphere models for several years. We find that the effect is 
less pronounced in J — if, where the difference is ~ 0.1. More extensive calculations of the 
TiO and H 2 line lists, while substantially improving the agreement with observed spectra 
at optical wavelengths, have both incr eased the discrepancy with the observed J — K color 
( lAllard et al.ll2000l ; iBaraffe et al.ll200ll ). Since TiO and H2O a re the dominant a bsorbers in 
the optical and near infrared, remaining limitations in the TiO (lAllard et al.ll2000l ). and more 



2 For brevity, we do not show the other near infrared CMDs that can be generated with the JHK band 
passes. The general appearance of the other diagrams is very similar to that of Mk vs J — K as far as the 
comparison between the modeled population and the data is concerned. 
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likely, the H2O (IPartridge fc Schwenkdll9971 ) line lists at high temperatures could reasonably 
be at the origin of this problem. On the other hand, dust opacity becomes significant in 
this regime and reddens the near infrared colors (Fig. [H]). At the T e g of late M dwarfs, the 
surface mass density of condensible material is modest, but the opacity of a tenuous cloud 
can be increased by decreasing the size of the grains. It is also possible that the onset of 
cloud format ion in late M dwar fs occurs at higher T e s than our models predict, as originally 
suggested by lTsuji et al.l (119961 ). The CMD suggests that the models may be underluminous 
rather than too blue (or a combination of both). The mass-ra dius relation for low mass stars 
is well-reproduced by evolution models down to ~ 0.096 M & (jChabrier &: Baraffell2000l ). but 
there is growing evidence that models of late L dwarfs are systematically underluminous for 
a given mass and age, as revealed by the first few dynam ical mass measurements below the 
stellar mass limit ([Ireland et al.ll2008l ; iDupuy et al.ll2008l ). For lat e T dwarfs, absolute mode l 
fluxes agree very well with observed spectral energy distributions (jSaumon et al.ll2006l 120071 ) 
so the theoretical radii of BDs older than > 1 Gyr and T c s < 900 K appear reliable. Finally, 
there are indications that the nearby field brown dwarfs may be a r e latively young po pulation 
(kllen et al.l I2OO5I : iMetchev fc Hillenbrand! I2OO6I : iLiu et al.l I2OO8I : ICruz et all l2008h . which 
would decrease the typical gravity compared to our simulation, increase the effect of clouds 
and move the population towards redder J — K (Fig. [8]). The discrepancy in the early L 
sequence could be the combination of a number of these effects. The input physics of the 
evolution model is fairly well understood and the comparisons between three independent 
calculations show that the results are fairly robust. On the other hand, the colors of model 
(e.g. J — K) are sensitive to the details of the input opacities of atmosphere models, which 
are known to be deficient in several ways. It appears more likely that the difference with 
the observed sequence arises from deficiencies in the atmosphere models rather than in the 
evolution. 



Small scale features in the synthetic CMDs are caused by the few cloudy atmosphere 
models that have converged to a poor solution, such as can be seen in Fig. [5J they have no 
physical significance. By construction, these CMDs fold in a distribution of gravities and 
T e ff based on reasonable assumptions that can be compared with the observed distribution. 
The concentration of the more massive objects along a narrow line in Fig. remains clearly 
visible in the synthetic CMD. This feature is not apparent in the data and it may be blurred 
by reasonable variations in metallicity (Fig. [TTj) The contribution of unresolved binaries to 
the blurring of the synthetic sequence is discussed in §4.2.1. The breadth of the late M and 
early L distribution is comparable to the scatter caused by variations in gravity, supporting 
the idea that the latter is i ndeed the "second parameter" of the sequence ( jTsuji fc Nakajima 
20031 : burrows et al.l 12006} ) . 



On the other hand, the widening of the observed distribution in J — K of the mid to 
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late L dwarfs cannot be explained solely by the range of gravity of the simulated population 
and appears to require variations in cloud properties. A combination of / se( j = 1 and 2 
populations appears sufficient to reproduce the observed range and the maximum J — K 
of late L dwarfs. This is confirmed by detailed spectral analysis of unusually red and blue 
J — K colors with our synthetic spectra, showing that L dwarfs sp an a range of / se d = 1 — 3 
(IBurgasser et al.l 120081 ; ICushing et al.l 120081 ; IStephens et al.l 120081 ) , with the redder dwarfs 
having lower values of / sec j. Variations in metallicity increasingly affect the J — K color of 
late L dwarfs and could also contribute to the widening of the observed L sequence (Fig. [TT]) . 
but probably not as much as variations in cloud properties. O bservations of binary L dwarfs 
such as the L4+L4 pair HD 130948BC JDupuv et al.l l2008h . will present a set of spectral 
standards of known mass, age, and metallicity that will allow such effects to be disentangled. 

Our simulations indicate that a third parameter, related to cloud opacity, is also nec- 
essary to model the distribution of the mid-L to mid-T dwarfs colors. The transition from 
L to T spectral types is thought to require a rapid change in the cloud properties over a 
small range of T e s. The most important free parameter in the lAckerman fc Marleyl (120011 ) 
cloud model is the sedimentation efficiency / se d- All cloud properties (median particle size, 
vertical structure, etc.) follow self-consistently given T eS , the gravity and the metallicity. 
Within this model, the transition from cloudy L dwarfs to the clear atmospheres of late T 
dwarfs can only be modeled by varying / se d. We find that the late L dwarf sequence requires 
fsed = 1 to 2, and that the transition dwarfs with 0.4 < J — K < 1.4 are fairly well matched 
with a / se d = 3 cloud (Fig. [TUb). Some parametrization of / sc d as a function of T cS (and 
perhaps gravity) could potentially match the CMD of L and T dwarfs. Admittedly, a color- 
magnitude diagram that involves only two band passes is a rather limited tool to study the 
transition. Detailed fits of the entire spectral energy distributions of ~ 20 LI to T6 dwarfs 
with these models fully support t he present hypothes i s as a trend of increa sing / sed for TO 
and later spectral types is found (ICushing et al.ll2008l ; IStephens et al.l 120081 ) . 



Finally, the J—K colors of late T dwarfs are best reproduced by our f se & = 4 model as our 
cloudless models are somewhat too blue. Adding a thin cloud layer makes the models slightly 
redder and brings them in better agreement with the data. Late T dwarfs are expected to 
be cloudless howeve r, as is indicated by detailed fits of their SED (jSaumon et al.ll2006l . 120071 ; 
Cushing et al.ll2008l ). The overestimated blue color of our cloudless models is most likely due 
to some inadequacy in the gas opacities, such as in the CH 4 line list or H 2 collision-induced 
absorption. In our models, the cloudless sequence is very narrow because the changes in Mk 
and J — K due to variations in gravity can be largely compensated by changes in T eS (Fig. [S]). 
Therefore, va riations in grav i ty can not explain the scatter in the colors of late T dwarfs, in 
contrast with Burrows et al. who find that gravities varying between log g = 4.5 and 

5.5 are adequate for the task. Our result is essentially independent of the gravity dependence 
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of the modeled colors. Late T dwarfs have 700 < T e ff < 1200 K. The (T e fj,logg) distribution 
of our a = 1 field population (Fig. [9]) shows that the majority (87%) of brown dwarfs in this 
T eff range have log g > 4.8, with a distribution that peaks sharply at log g = 5.3. Objects 
with logg ~ 4.5 nearly all have cooled well below the T e g of any object detected so far. 
Unless the sample of BDs with parallaxes contains a large fraction of objects that are much 
younger than 5 Gyr (the average age in this simulation), there are too few low gravity objects 
in the late T dwarf T e fj range to account for the observed spread of the sequence , even when 
considering the different JHK colors of the cloudless sequence of iBurrows et all (120061 ). It is 
conceivable that some late T dwarfs may not have completely cloudless atmospheres, which 
would increase the variation in near-infrared colors of late T dwarfs. An example of such 
an objec t is the T5.5 dwarf S PSS Jll 1009.99+011613.0 whose SED favors a model with 
iscd = 4 ([Stephens et al.ll2008l ). The dispersion caused by metallicity variations is discussed 
in 84.2.1. 



4.2. A hybrid sequence of models 



To further illustrate the potential of population synthesis for the study of brown dwarfs, 
we have developed a simple model of the L/T transition. Based on the above discussion of 
the near-infrared CMD, we consider an evolution sequence where the cloud sedimentation 
parameter / sed increases as a function of decreasing T eff . We linearly interpolate the surface 
boundary condition in T e g- (Fig. [I]) between cloudless models at 1200 K and cloudy models 
(/sed — 2) at 1400 K for each gravity. For simplicity, we assume that the characteristics of the 
transition depend only on T e g. The role of gravity in the physics of the transition, if any, will 
have to be determined empirically. We use the magnitudes of / se d = 4 atmosphere models 
for T e fj < 1200 K and of / sec j = 1 models above 1400 K. The slight inconsistencies between 
the / S ed values used in the evolution sequences and for computing the magnitudes are modest 
and can be ignored here (see §2.5). The colors are linearly interpolated in T e g at constant 
gravity between those two regimes. Our linear interpolation between cloudy and cloudless 
models implies a specific but arbitrary choice of the cloud evolution with T e g a cross the 
L/T t ransition. This approach to modeling the L/T transition is simi lar to that of IChabrier 



( 2002 ) who inte rpolated magnitudes in T e ff between the DUSTY00 (IChabrier et al. 



and the COND (IChabrier fc Baraffdl2000l ) evolution sequences and to lBurgasser et al. 



2000a[) 

faoj ' 



who used a slightly more elaborate interpolation between cloudy and cloudless magnitudes 
to approximate the effect of cloud clearing at constant T e g. Here we have recomputed the 
cooling with the hybrid bo undary conditio n wh ich results in fea t ures in the evolution that 
could not occur within the IChabrier! (120021 ) and iBurgasser et all ( 120021 ) schemes. 
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The (T e s,logg) distribution for this hybrid evolution sequence, under the same IMF 
and SFR assumptions as in Fig. [9] is shown in Fig. [12j Two new features are apparent. 
The most apparent is a bump in the upper envelope of the distribution, where the maximum 
gravity rises at the transition from the lower value of cloudy models to the higher value 
of cloudless models, as discussed in §3.1. This effect is modest (~ 0.1 in logg) but is a 
direct consequence of the change in surface boundary condition as the atmosphere becomes 
free of clouds. In reality, this bump may be narrower or more spread out than is modeled 
here. The other feature is an increased density of brown dwarfs in the 1400-1200 K range 
for log g > 4.8. The relatively high entropy interior (equivalently, higher Tio, see Fig. [I]) 
of a 1400 K cloudy brown dwarf must release much more of its heat content to become 
a 1200 K cloudless brown dwarf than would normally be required for cooling to the same 
effective temperature without the transition in cloud properties. This increases the cooling 
time scale, i.e. the brown dwarf cooling slows between 1400 K and 1200 K and causes a 
moderate pile up in the T e g range of the transition. The effect can be clearly seen in the 
T e fj distribution of the hybrid disk sequence where an excess of brown dwarfs of a factor of 
2.1 occurs in the transition region compared to either the cloudy or the cloudless evolution 
sequences (Fig. [T3]). For the same reason, an excess also occurs for the cloudy sequence 
relative to the cloudless case for 650 < T e s < 1000 K as the cloud sinks in the atmosphere 
to levels where it no longer affects the surface boundary condition (Fig. [Q. This is a more 
gradual effect than our modeled L /T transition and the excess brown dwarf density is spread 
over a wider range of temperatures. Qualitatively, this effect is an inevitable consequence of 
the disappearance of clouds at the L/T transition and is not dependent on any particular 
physical mechanism or model of the "cloud collapse" responsible for the transition. 



This excess of brown dwarfs in the transition is in apparent contradiction with lBurgasser 



(120071 ) who found a minimum in the space density of transition dwarfs as a function of spectral 
type. There are two reasons for this. Our simulation of cloudy and cloudless evolution 
indicates a m onotonous i ncrea se in the T e g distribution towards lower temperatures. As 



pointed out by lBurgasserl (120071 ). once converted to a distribution in terms of spectral type (or 
absolute magnitude) the distribution would show a minimum because the relation between 
Tpff and spectral type is highly non-linear and shows a plateau in the L/T transition region 



(jGolimowski et al.ll2004l ). The other effect is the slowing of the evolution of brown dwarfs 
in our hybrid evolution sequence as the brown dwarf converts from a cloudy to a cloudless 
atmosphere. To the best of our knowledge, this effect is not included in any of the earlier 
population syntheses which are based on a single evolution sequence (cloudy or cloudless). 

Also shown on Fig. dH] are (T e fj,logg) determinations for various field and ultracool 
dwarfs in binaries. Given the small number of constrained objects and the current lack of 
detections at very low T e g, the agreement between the modeled and observed populations is 
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reasonable. As spectral fitting becomes more precise, particularly for cloudy L dwarfs and 
for binary brown dwarfs, and with the creation of a larger, volume-limited sample of BDs 
with parallaxes, such comparisons will test our understanding of ultracool dwarf formation 
and evolution. 

The corresponding CMD is shown in Fig. [TH The general agreement with the data is 
quite good as we have chosen the parameters of the hybrid model sequence for that purpose. 
Figure [His to be compared to Fig. [TOh, for T cS > 1400 K and Fig. 0TOH for T cff < 1200 K. The 
region of interest is the transition region where the synthetic population shows a well-defined 
lower envelope whose shape and location depend on the details of the transition process. 
The existence of this lower boundary simply reflects the high gravity envelope of Fig. [12] 
There is increased scatter in the transition region that arises from the distribution in gravity 
and cloud sedimentation parameter. The larger density of brown dwarfs at T c g ~ 1300 K 
seen in Fig. [12] is also visible in the CMD, centered at J — K ~ 1. The clump of five 
early T dwarfs is superposed on this density enhancement in the synthetic CMD but this 
is not a statistically significant result. This subtle feature may be discernible in a much 
larger data set, however. Adaptive optics imaging of early T dwarfs has res olved many into 



earlier and later comp onents, effectively depopulating the transition region (ILiu et all 12006 



Burgasser et al.ll2006l ). Nonetheless, the cooling of brown dwarfs is a continuous process and 
the full range of T e g must be populated regardless of the distribution as a function of spectral 
type. Our synthetic CMD shows that there should be no significant dearth of objects in the 
transition region (J — K ~ 0.2 - 1.5). The current paucity of known transition objects must 
be a selection effect. The above remarks also apply to the Mj vs J — K CMD where we 
obtain an equally good agreement with the data, including the observed brightening in the J 
band. The Mh vs J — H CMD shows the same general features but with a larger dispersion 
below the transition. While the transition can be modeled with success by varying / se d in our 
atmosphere models, we do not propose that this is the physical mechanism of the transition. 
It is merely a parametric description of the transition using a specific (and simple) cloud 
model. 



4-2.1. Synthetic CMD with different assumptions 

In this section, we consider different assumptions for our population synthesis with the 
hybrid evolution sequence to estimate the importance of variations in the SFR, the IMF, 
and the role of unresolved binaries. These experiments are summarized in Fig. [15] where our 
fiducial simulation of Fig. [12] is reproduced in panel a for convenience. Since ther e is evidence 



that the local BD population has a younger average age than solar-type stars (lAUen et al. 
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2005 : iMetchev k Hillenbrand! 120061 : IZapatero Osorio et all 120071 : iLiu et al.ll2008l : ICruz et al. 
20081 ). we consider a SFR where all BDs were formed at a constant rate in the past 5 Gyr 
(Fig. IT5"b). This results in a larger fraction of low-gravity objects, widening the spread of 
the sequence all the way to J — K ~ 0. The sharp blue edge of the sequence remains well- 
defined above Mk = 12 but the wider sequence reproduces the data fairly well. If indeed, 
the local population of BDs turns out to be younger than our fiducial SFR would imply, this 
would reduce the contribution of cloud properties as the "second parameter" favoring the age 
(gravity). In this case, the widening of the late L sequence is caused primarily by the turnover 
in colors at the onset of the transition. While the observational evidence is consistent with a 



roughly flat SFR over the age of of the galaxy (jAUen et al.ll2005l ) we consider an expon entially 
decreasing SFR with a characteristic time scale of r = 5 Gyr (jMiller &: Scald Il979l ). This 
results in a larger fraction of old BDs and a very narrow distribution in the CMD (Fig. 
[T5b ) that bears little resemblance with the data. In panel d, we show the effect of using 



the log-normal IMF of IChabrier et al.l (120051 ) which greatly reduces the number of low mass 
objects compared to our fiducial power law IMF with a = 1. This reduces the number of 
low T e g and low-gravity objects and the sequence becomes somewhat narrower for Mk > 13. 

The i mpor tanc e of binaries in bro wn dwarf population synthesis has been pointed out by 



Chabrier 


(2002 


) and 


Bureasser 


(2007) 



resolution imaging surveys aimed at establishing the statistics and properties of binary brown 
dwarfs. While the number of resolved systems remains small (~ 30), the observed binary 
fraction is Estimating the contribu tion of unresolv ed systems, the binary fraction 



is somehwere in the range e& ~ 0.1 — 0.6 (jBurgasser 2007 ). The mass ratio q = M 2 /M± 
distribution stron gly favors equal mass binaries and can be approximated by f(q) ~ q A 
( iBurgasserl (120071 ) and references therein). For the purpose of illustrating the role of binaries 
in the synthetic CMD we take a fraction = 0.3 of the simulated single star population and 
assign each of them a companion of the same age and with a mass drawn from the above mass 
ratio distribution. The magnitudes of the components are combined to get the magnitudes of 
the (assumed) unresolved binary. The result is shown in Fig. [T5e . where unresolved binaries 
are shown in light blue and single dwarfs in green. Because the mass ratio distribution is 
sharply peaked at q — 1, the binaries systems form a sequence that is approximately that of 
single dwarfs shifted upward by 0.75 magnitude. The sharp boundaries of the sequence (on 
the blue side for Mk < 12 and on the faint side of the transition region) remain well defined 
even though the total dispersion has increased significantly. Unless the binary mass fraction 
is well under our choice of e& = 0.3 and the mass ratio distribution is much flatter than g 4 , 
unresolved binaries contribute significantly to the breadth of the dwarf sequence at all T eS . 



Finally, we consider the effects of variations in metallicity on the CMD. The large grids 
of cloudy atmosphere models with non-solar metallicity required for such a calculation are 
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not yet available. We can presently explore the effect of metallicity on cloudless atmosphere 
models with [M/H]=-0.3, and +0.3. These are appropriate for T dwarfs beyond the 
L/T transition. Separate evolution sequences were computed for each of those three values 
of [M/H]. In addition to the input mass and age distributions, the sy nthetic population 



is dr awn from the observed distribution in metallicity of 253 M dwarfs (ICasagrande et al. 



20081 1 spanning -0.9 < [M/H] < 0.3 with and average of -0.20 ± 0.22. Our model grids 



encompasses 73% of the stars in this distribution so extrapolations of the evolution and 
colors to metallicities below —0.3 play a small role. For this calculation, we adopt a power 
law IMF with a = 1 and a constant SFR for ages of 0-10 Gyr, as in Fig. [TUb . As we found for 
cloudless models with solar metallicity, the distribution is systematically too blue in J — K 
by about 0.4 magnitudes, but the slope and width of the distribution is quite similar to 
that of the observed late T dwarf sequence. Thus metallicity is likely to be a significant 
contributor to the "second parameter" for the late T dwarfs. 



4.3. Color-Magnitude diagrams for clusters 



Brown dwarfs in galactic clusters offer the advantages of a uniform population with a 
known distance, metallicity and a very narrow age distribution compared to field objects. 
On the other hand, brown dwarfs in clusters are typically much more distant than typical 
field brown dwarfs and there are very few cluster L and T dwarfs known because they are 
extremely faint. They could als o be rare i f the I MF decreases towards lower masses, as with 
the log-normal distribution of IChabrierl (120021 ). Technological progress will undoubtedly 
help mitigate these limitations and we anticipate the emergence of new diagnostics of the 
evolution of brown dwarfs from the study of galactic clusters. We have generated synthetic 
populations of 1000 objects with masses between 0.006 and 0.1 M & and an IMF index of 
a = 0.6 and for ages of 10, 50, 100, 200 and 5 00 Myr, with our hyb r id evolution sequence . 
This value of a is typical for galactic c lusters iMoraux et al.l (120031 ) ; iBouvier et al.l (120031 ) ; 
Barrado y Navascues et al.l (120021 . 12004 ) , although a re cent study suggests th at the IMF of 
young clusters is decreasing in the substellar regime ([Andersen et al.l 120081 ) . The SFR is 
assumed constant for ±1 Myr around the age of each cluster, except for the 10 Myr case 
where we used a dispersion of ±0.5 Myr. The resulting (T e ff, \ogg) distributions, which 
follow isochrones of the hybrid evolution sequence, are shown in Fig. [161 The phase of 
deuterium burning is particularly prominent at these relatively young ages and occurs within 
the modeled transition for ages of 50-200 Myr. Inflexions in the isochrones at the T e g- end 
points of the modeled L/T transition are also apparent. 



All five synthetic clusters are shown in the Mk vs J — K CMD (Fig. [17]) . where objects 
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with T fj > 2400 K have been removed as they fall outside of our surface boundary condition 
in a regime where extrapolation is unreliable. This figure is to be compared with the CMD for 
the field population (Fig. [Tij) . In relatively young clusters (10 - 50Myr), the L/T transition 
region occurs at 1.5 to 2 magnitudes brighter in than for field brown dwarfs, on account 
of their lower gravity. The J — K of the reddest brown dwarfs in a cluster increases steadily 
with the cluster age, under our assumption that the T e g of the onset of the transition is 
independent of gravity. 

The most unusual (and to our knowledge, not previously predicted) feature of the cluster 
CMD is a backtracking loop along each isochrone except the oldest (500 Myr). It is not 
plotted on the youngest sequence (10 Myr) because it occurs at T eS ~ 2600 K, above the 
limit of our evolution tracks. This feature arises because at ages greater than about 10 6 yr, 
brown dwarfs in the deuterium burning phase (which occurs at lower T e g and larger Mk 
for older clusters (Fig. [TBI) ) are more luminous than both less massive and somewhat more 
massive objects. It is a simple consequence of the multi-valued gravity (i.e. luminosity) for a 
given T e fj (i.e. color) around the phase of deuterium burning (Fig. [16]). This feature should 
be most distinctive and detectable in clusters with ages of the order of ~ 50 to 10 Myr. This 



featur e is also visible as a bump in the cluster LF, as seen in the simulations of lAUen et al. 



(12003( 1 



4.3.1. Pleiades 



We compare our simple model of the L/T transition with the results of deep brown 
dwarf surve ys in the Pleiades clu ster. The relevant parame ters of the clu s ter ar e its distance 
(134 ± 3 pc, IPercival et aD j2005h. age (~ 100 - 120 Myr, iMartm et all (jl998h ). metallicity 



f[Fe/Hl= -0.034 ± . 024, 



Boesgaard fc Friell (119900 ). and extinction (E(V — I) 



0.06, 

Stauffer fc Hartmannl (119871 )). Deep surveys have established that its low mass stellar and 
substellar population can be reproduce d by a power law IMF with an index of a ~ 0.6 



jMoraux et al.l I2OO3I : iBihain et alll2006h . We generate a synthetic brown dwarf population 
with our hybrid evolution sequence QM/H] = 0), with a = 0.6 and stellar ages uniformly 
distrib uted between 105 an d 11 5Myrn The synthetic cluster sequence is compared with the 



data of IBihain et al.l (120061 ) and lCasewell et al.l (120071 ) in Figs. [181 and [T9l r espec tively, where 
we assume that E(J — K) = A K = and d = 134 pc. The lCasewell et al.l (120071 1 survey goes 
> 1 magnitude fainter and shows a tighter L dwarf sequence. Both figures show that the 



3 This probably overestimates the age spread of the cluster but this hardly matters here, as shown by the 
narrow sequence of the synthesized population in Fig. 1181 



-23- 



modeled L dwarf sequence is too bl ue by ~ 0.3 magnitude , too faint by ~ 0.5 magnitude, or 



a combination of both effects. The ICasewell et al.l (120071 ) data (Fig. [19]) goes deep enough 
to reveal two objects that appear to be early T dwarfs in L/T transition for the cluster, but 
they are a good magnitude fainter than the modeled population. Given that this hybrid 
evolution/color model agrees fairly well with the local field population of brown dwarfs (Fig. 
[141) . the discrepancy shown in Figs. [18] and [19] is somewhat surprising. 



The part of the modeled sequence that corresponds to the M and L spectral types 
[Mr < 12) is ~ 0.4 magnitude too blue in J — K. A younger age for the sequence would 
improve the agreement but a cluster age of 30 — 40 Myr is required to make it overlap the 
data - an implausible explanation. Figure [20] shows that this departure from the data is not 
much worse than we found for the disk population (§4.1) and in all probability is caused by 
the sa me model deficiency. T he two dimmest objects that appear to be in the L/T transition 
region (ICasewell et al.l 120071 ) are almost one magnitude fainter than predicted by our simple 
transition model. These T dwarf candidates could be matched by lowering the T e g regime 
of the transition, which is set to 1200 — 1400 K in our hybrid sequence to match the CMD 
of field brown dwarfs. These two Pleiads suggest the possibility that the onset of the L/T 
transition occurs at lower T eS in lower gravity objects. This hypothesis could be tested with 
a detailed study of a larger sample of dim Pleiades brown dwarfs in the range of K ~ 17—19, 
for example. It is supported by the low T e gf fou nd for two transition brown dw arf companions 
of young main seque nce stars, HD 203030 B ( iMetchev fc Hillenbrand] 120061 ) and HN Peg B 
(ILuhman et al.ll2007l ). This may be an important clue as to the physical trigger of the "cloud 
collapse" that characterizes the transition. On the other hand, the association of these two 
brown dwarfs with the Pleiades is uncertain due to rather large error bars in their parallaxes. 
Interpretation of the mismatch should be regarded with caution until membership can be 
confirmed. Follow-up spectroscopy should unambiguously reveal their nature. Finally, the 
phase of deuterium burning causes a well-defined feature in the transition region of the near- 
infrared CMDs that is potentially observable in this cluster if more members with K ~ 16.5 
- 18 can be identified. In the field, this feature is not discernible because all the primordial 
deuterium has been consumed in all but the very lowest mass objects for ages above 1 Gyr 
(Figs. [2]and[l]). 

Figure [2IH combines both the field and Pleiades data as well as the corresponding modeled 
sequences. The Pleiades M and L sequence objects are brighter than the those in the field due 
do their larger radii at a much younger age. Above the transition, the modeled sequences 
are both systematically too blue (or too faint, but that is less likely) but show the same 
behavior as the data with the same approximate variation between the two populations. In 
both cases, the maximum J — K reached as well as the slightly redder value for the older field 
population is well reproduced. Interestingly, the two faintest Pleiades candidates appear to 
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be in the L/T transition fall within the field dwarf transition region. 



5. Conclusions 

Our calculation of the evolution of very low mass stars, brown dwarfs and planetary 
mass objects produces models that are quantitatively in very good agreement with pub- 
lished calculations. A detailed comparison shows some systematic differences that can be 
attributed to conductive energy transport, different choices of composition for the interior, 
and for the initial state but the primary source of discrepancy is the surface boundary con- 
dition. The cloudless and cloudy model atmospheres from which we extract the surface 
boundary condition have been validated with extensive comparisons with spectroscopic and 
photometric data. Thus, our boundary condition is quite realistic and we do not expect that 
the foreseeable improvements in atmosphere models will have much effect on our modeled 
evolution of brown dwarfs. By using the boundary condition from our atmosphere models, 
we can compute self-consistently the evolution, absolute fluxes, absolute magnitudes, and 
colo rs for a variety of clo ud properties, metallicities and eventually including vertical mix- 



ing (jSaumon et al.l 120061 ) . Only the evolution across the L/T transition region cannot be 



modeled in this self-consistent fashion because a transition cloud model is not yet available. 

We have developed a simple model for the cooling and color evolution of brown dwarfs 
across the L/T transition. This hybrid model predicts an excess of brown dwarfs in the 
T e ff range of the transition by about a factor of 2 compared to purely cloudy or cloudless 
evolution. We have applied this hybrid evolution model, combined with the near-infrared 
magnitudes predicted by our atmosphere models, to generate synthetic CMDs that can be 
compared with samples of brown dwarfs in the field and in galactic clusters. Our primary 
focus is the "second parameter" responsible for the dispersion about the brown dwarf se- 
quence in the CMD. Population synthesis is a potentially powerful tool but the results can 
be quite sensitive to the input assumptions for the IMF, SFR, metallicity, and the binary 
frequency and mass ratio distribution. Our knowledge of the brown dwarf population in the 
solar neighborhood cannot yet provide these inputs a priori. The relatively small number of 
brown dwarfs with known parallax and the heterogeneous nature of the sample imply that 
the observed distribution of field substellar dwarfs in the CMD is not yet fully characterized. 
Nevertheless, both observations and models have reached a stage where general trends in 
the "second parameter" along the L-T spectral sequence can be interpreted. 

We find that for our fiducial assumptions (power law IMF with a — 1, constant SFR over 
the past 10 Gyr, single brown dwarfs only and [M/H]=0), the hybrid sequence reproduces the 
overall sequence from late M through late T dwarfs rather well, but not the dispersion along 
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the sequence. Based on the near-infrared CMDs, we find that the L/T transition occurs 
between T e g ~ 1400 — 1200 K in field brown dwarfs, in agreement with previous estimates. 
While a transition over such a narrow range of T e g appears to be "fast" considering the 
rather dramatic change in JHK colors across the transition, these values of T e g correspond 
to ages of 2 and 4 Gyr for a 0.06 M & brown dwarf. The duration of the transition decreases 
rapidly with mass however, lasting only 0.15 Gyr for a 0.03 M® brown dwarf. 

Better agreement can be obtained from late M to late L spectral types if the population 
is younger, such as with a constant SFR that started only 5 Gyr ago, by including binaries, 
or assuming that there is a wider range of cloud properties for later L spectral types. For a 
fixed metallicity, all simulations predict that the distribution of brown dwarfs in the CMD 
will have a sharp edge formed by old brown dwarfs of all masses (> 3 Gyr). This edge is 
on the blue side of the distribution for Mr < 12.5 and to the red side after the J — K 
color of the sequence turns over, corresponding to Mk > 13. This feature is not visible in 
the data, however, most likely because it is blurred by variations in metallicity within the 
sample. We are not able to include metallicity variations in simulations of cloudy brown 
dwarfs, except in a very approximate way (Fig. [TT1) . We find that it could be a significant 
contributor to the second parameter. Detailed spectral analysis of brown dw arfs with unusual 



J — K colors for their spe c tral types have more extreme cloud parameters (IBurgasser et al. 



20081 ; ICushing et al.l 120081 ; IStephens et al.ll2008l ). which strongly suggest that cloudiness is 



the second parameter. Spectral analysis with models of non-solar metallicities have barely 
begun, however, so it would be premature to attribute the dispersion along the L sequence 
entirely to cloud characteristics. On the other hand, a simulation of cloudless models with 
an empirical metallicity distribution shows a good match to the dispersion of the late T 
dwarfs. For those coolest dwarfs, we find that the second parameter is a combination of 
metallicity variations (which dominate) and binaries. Gravity is not important as the with 
of the distribution is not affected by the choice of SFR or IMF. To summarize, we find 
that there is no single second parameter that accounts for the dispersion of brown dwarfs 
around the sequence seen in the CMD. A young age distribution, a range of metallicities and 
cloud properties as well as binaries all contribute to the dispersion. The challenge will be to 
untangle their contributions. 

The hybrid model fares somewhat worse when compared to the much younger brown 
dwarf population of the Pleiades. If the two faintest Pleiads reported are indeed T dwarf 
members of the cluster, then they provide strong evidence that the L/T transition occurs 
at lower T e g in lower gravity objects (i.e. younger or less massive). Finally, isochrones in 
CMDs clearly reveal the phase of deuterium burning at young ages, a feature that should be 
observable in young clusters with ages between of ~ 50 — 100 Myr. 
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At this time, the modest size of the sample of L and T dwarfs with known parallax and 
the lingering problems in modeling the atmospheres of cloudless and cloudy brown dwarfs 
restrict how much we can learn from the study of CMDs. Model limitations will eventually be 
overcome as new moelcular line lists are being developed for key molecules and cloud models 
become more sophisticated. We anticipate a rich harvest of brown dwarf parallaxes from the 
volume limited solar neighborhood census component of the Panoramic Survey Telescope 
& Rapid Response System (Pan-STARRS) @. Color-magnitude diagrams are a potentially 
powerful tool for the study of brown dwarf evolution and of the L/T transition. Statistical 
comparisons with synthetic populations in two-dimensional parameter space will become an 
important complement to the detailed studies of the spectra of individual transition objects 
and of brown dwarf binaries. 
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and J. R. Stauffer for useful discussions. 



A. Typographical errors in Saumon, Chabrier & Van Horn (1995) 

Typographical errors have been found in the expressions for the thermodynamics of 
hydrogen and helium mixtures given in Saumon, Chabrier & Van Horn (1995, hereafter 
SCVH). The correct expressions are: 



C _ f-i v\^ H oH \^^ He cHe , $naix 9 log S m j, 

S t-(1-Y)^S t +Y—S t +—^-jr 



(45) 



(46) 



3(1 — Vh 2 — Vh) 
where the equation numbers are those of SCVH. 



(56) 



The last correction affects the calculation of the contribution of the electrons to the ideal 
entropy of mixing. In practice, this matters only when the gas is nearly fully ionized, where 



^http : //pan-starrs . if a.hawaii . edu/project/reviews/PreCoDR/documents/scienceproposals/ sol .pdf 
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the ideal entropy of mixing S m ^ (Eq. 53) is only a rough approximation of the actual entropy 
of mixing, except in the high temperature, low density limit where it is nearly exact. This 
correction should not be of much concern to applications of the mixed H/He EOS to stellar 
and planetary interiors where only a very small fraction of the mass is not fully ionized. 
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Fig. 1. — Surface boundary condition, expressed as T 10 (T eff , g) for solar metallicity atmo- 
spheres without clouds (dotted curves) and with clouds, using a sedimentation parameter 
f se d = 2 (solid curves). Each curves correspond to a different surface gravity with logg = 4, 
4.5, 5, and 5.5, from top to bottom, respectively. Larger values of T 10 imply a higher entropy 
for the interior. The surface boundary condition is interpolated below T cS = 500 K to (0, 0). 
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Fig. 2. — Evolution of the luminosity of brown dwarfs for masses of 0.005, 0.01, and up to 
0.08 M in steps of 0.01, from bottom to top, respectively. This sequence of models uses 
cloudless model atmospheres for the surface boundary condition. Solid dots show the models 
of the Lyon group based on the COND03 atmosphere models for the same masses (Baraffe 
et al. 2003). 
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Fig. 3. — Same as Fig. [2] but with a cloudy surface boundary condition (/ se d = 2). Solid dots 
show the models of the Lyon group based on the DUSTYOO atmosphere models (Chabrier 
et al. 2000; Baraffe et al. 2001). 
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Fig. 4. — Evolution of brown dwarfs in T e ff and gravity for the sequence based on cloudless 
atmospheres. The evolution proceeds from right to left along the heavy black lines, which 
are labeled with the mass in M . Isochrones are shown by the blue dotted lines: (from right 
to left) 0.01, 0.02, 0.04, 0.1, 0.2, 0.4, 1, 2, 4, and lOGyr. The nearly vertical red lines are 
curves of constant luminosity: (from right to left): logL/L = —3, to —6.5 in steps of 0.5. 
Curves of constant radius are shown in green: (from top to bottom) 0.08 to 0.13 R & in steps 
of 0.01. The phase of deuterium burning is revealed by the kink in the isochrones for objects 
with masses between 0.01 and 0.02 M & . [See the electronic edition of the Journal for a color 
version of this figure.} 
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Fig. 5. — Same as Fig. H]but with cloudy atmospheres models for the surface boundary con- 
dition (f se d = 2). The cloudless cooling tracks from Fig. [Hare shown in gray for comparison. 
[See the electronic edition of the Journal for a color version of this figure.} 
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Fig. 6. — Evolution of the luminosity of brown dwarfs This shows the same clo udless sequence 
as in Fig. [2] (solid lines), compared with the cloudless cooling sequences of iBurrows et al. 
Jl997h (dotted lines). 
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Fig. 7. — Col or-magnitude d i agram (M KO system) f o r field brown dwarfs, showing the 
photometry of iLeggett et al.l (120021 ) and iKnapp et al.l (12004 ) with M dwarfs in black, L 
dwarfs in red and T dwarfs in blue. All known binaries h ave been removed from the sample 



except those with resolved MKO photometry: e Indi B ( McCaughrean 



102109.69-030420.1 and SDSS ,1042348.57-0 41403.5 (IBurgasser et al 



..Liu &: Leggettll2005l). The para l laxes are from iPerryman et al.l ( 
Tinney et al.~(j2003 ). Vrba et al. ( 2004 ). and various references in 



199 



71. 



et al. 120041 ) . SDSS 



2006), and Kelu-1 



Dahn et al 



Leggett et al.l (120021 ) . [See 



20021 ). 



the electronic edition of the Journal for a color version of this figure] 
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Fig. 8. — Color-magnitude diagram (MKO system) for atmosphere models at three gravities 
(log g = 4.5, 5 and 5.5, from right to left at the bright end) and for cloud condensation 
parameters / se d = 1 (red), 2 (blue), 3 (green), 4 (cyan), and cloudless models (black). The 
effective temperature ranges from 2400 K (top) to 500 K and squares along each curve indicate 
T e Q = 2000, 1500, and 1000 K (from top to bottom). All models shown have solar metallicity. 
Some of the curves are not smooth due to numerical difficulties with a few cloudy cases. [See 
the electronic edition of the Journal for a color version of this figure.} 
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Fig. 9. — Distribution in (T e s,\ogg) of a simulated local field population of brown dwarfs 
with power law IMF with index a = —1, masses between 0.006 and 0.1 M , and a uniform 
age distribution between and 10 Gyr. The mapping from (M,age) to (T eff , logg) is based 
on the cloudy evolution sequence with / scd = 2 and [M/H]=0. 
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Fig. 10. — Color-magnitude diagram of a synthetic field brown population of solar metallicity 
based on the cloudy evolution sequence (f se d = 2) and cloudy model magnitudes with / sed = 
1—4 (panels a-d) and for the cloudless sequence (panel e). The synthetic population is taken 
from Fig. [9j The sources of data are the same as in Fig. [7J The small green squares show 
the synthetic population. All magnitudes are in the MKO system. [See the electronic edition 
of the Journal for a color version of this figure.] 
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Fig. 11. — Colors of cloudy models with / se d = 2 and logg = 5 and two different metallicities 
(solid lines). The curve on the right has [M/H]=+0.5 and extends from T e g = 2400 to 
1300 K. The curve on the left, which also appears on Fig. [HJ has [M/H]=0 and extends from 
T e fj = 2400 to 500 K [See the electronic edition of the Journal for a color version of this 
figure.] 
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Fig. 12. — Distribution in (T e g,logg) of a simulated local field population of brown dwarfs 
with power law IMF index a — 1, masses between 0.006 and 0.1 M Q , and a uniform age distri- 
bution between and 10 Gyr. The mapping from (M,age) to (T e g, log (?) is based on a hybrid 
cloudy/cloudless evolution sequence to mimic the L/T transition (see text). The vertical lines 
indicate the T e g range of the modeled transition. Superimposed on the synthetic distribution 
are the (T e fi,\o K q) of field dwarfs determined by fitting their spectral energy distribution 
(ICushing et al.l (120081 ). blue). The red boxes show the parameters of well-studied dwarfs 
that are constrained by L and age. From left t o righ t, they are 2 MASS J041519 5— 093506 



(ISaumon et all l2007h. HP 3651B flLiu et all [20071), Gl 570D flSaumon et al 



200 



Ind Bb (IMcCaughrean et al. I2004T). 2 MASS J1 2171110- 03111131 flSaumon et all 120071) 



2MASS J09373487+ 2931409 flGeballe et al.ll2008[). Gl 229 B (jSaumon et al.lboooh . HN Peg B 



(ILeggett et al.ll2008h . e Ind Ba (IMcCaughrean et al.ll2004), and 2 MASS J05591914-1404488 



(logL/L = —4.606 ±0.04 obtained with the method of iGeballe et al.l (120011 ) assuming it is 
a single brown dwarf). For objects without age constrains available, 1 — 10 Gyr is assumed. 
[See the electronic edition of the Journal for a color version of this figure.} 
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Fig. 13. — Distributions in T cS of the synthetic disk population for the hybrid sequence 
(black), the cloudy sequence (red) and the cloudless sequence (blue). The distributions are 
normalized to the size of the sample. The sharp cutoff at 200 K reflects the mass cutoff of our 
simulations at 0.006 M Q . By construction, the hybrid sequence recovers the cloudy sequence 
above T eS = 1400 K and the cloudless sequence below 1200 K. In the hybrid sequence, the 
cooling time scale increases in the transition region (1200-1400 K) causing a pile up of objects. 
A similar effect occurs in the cloudy sequence when the cloud deck becomes so deep that it 
no longer affects the evolution. The result is an excess of objects in the 600-1000 K range. 
[See the electronic edition of the Journal for a color version of this figure.] 
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Fig. 14. — Synthetic color-magnitude diagram (MKO system) of a field brown dwarf pop- 
ulation based on the hybrid evolution sequence and colors (see text). See Fig. [TU] for the 
legend. [See the electronic edition of the Journal for a color version of this figure] 



-47- 




-10 1 2-10 1 2-10 1 2 
J-K J-K J-K 



Fig. 15. — CMD of the simulated field population resulting from variations in the input 
assumptions, a) fiducial sequence based on the hybrid evolution sequence with a power law 
IMF (a = 1), a flat SFR over 0-10 Gyr, and single stars only (same as Fig. [T4l) b) same as 
in a) but the SFR is truncated for ages > 5 Gyr. c) same as a) but with a n exponentially 



decrea sing SFR (r = 5 Gyr). d) same as a) but with the log-normal IMF of IChabrier et al. 



( 120051 ) . e) same as a) but with unresolved binaries included (light blue dots). The binary 
mass fraction is e& = 0.3 and the distribution of mass ratios is f(q) ~ q A . f) cloudless models 
with the metallicity distribution of M dwarfs (ICasagrande et al.ll2008l ). a constant SFR over 
0-10 Gyr and a power law IMF with a — 1. [See the electronic edition of the Journal for a 
color version of this figure.] 
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Fig. 16. — Distribution in (T eff ,log g) of simulated cluster populations of brown dwarfs with 
power law IMF index a = 0.6, masses between 0.006 and 0.1 M , and ages of 10 (black), 50 
(red), 100 (blue), 200 (green) and 500 Myr (magenta). The mapping from (M,age) to (T e g, 
log p) is based on a hybrid cloudy /cloudless evolution sequence to mimic the L/T transition. 
The vertical lines indicate the T e fj range of the modeled transition. [See the electronic edition 
of the Journal for a color version of this figure.] 
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Fig. 17. — Synthetic color-magnitude diagram (MKO system) of cluster brown dwarf pop- 
ulation shown in Fig. [161 Each sequence corresponds to a different age: 10 Myr (black), 
50 Myr (red), 100 Myr (blue), 200 Myr (green), and 500 Myr (magenta), from top to bottom, 
respectively. [See the electronic edition of the Journal for a color version of this figure] 
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Fig. 18. — Color-magnitude diagram of the Pleiades. A brown dwarf population based on 
the hybrid evolution sequence and colors (see text) with an IMF index of a = 0.6 and an 
age of 110 Myr with a ±5Myr dispersion is shown in green. The black squares show the 
2MASS photometry reported in lBihain et al.l (120061 ) for brown dwarf candidates with cluster 
membership confirmed by proper motion measurements. [See the electronic edition of the 
Journal for a color version of this figure] 
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Fig. 19. — Color-magnitude diagram of the Pleiades. A brown dwarf population based on 
the hybrid evolution sequence and colors (see text) with an IMF index of a = 0.6 and an age 
of 110±5 Myr i s show n in green. The black squares show the MKO photometry reported in 
Casewell et al.l (120071 ). [See the electronic edition of the Journal for a color version of this 
figure.} 
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Fig. 20. — Color-magnitud e diagram (MKO system) showing brown dwarfs in the Pleiades 
(red, ICasewell et al.l (120071 )). and in the field (blue, see the caption of Fig. [7] for references). 
The corresponding hybrid model sequences are shown by small dots in lighter colors. [See 
the electronic edition of the Journal for a color version of this figure.} 



